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RETURN-CURRENT-DRIVEN  INSTABILITIES  OF  PROPAGATING  ELECTRON  BEAMS 

I.  Introduction 

In  recent  years,  there  have  been  many  investigations *~^of  the  resistive 

instabilities  of  a  self-pinched  relativistic  electron  bean  propagating  through  a 

collisional  plasma.  These  instabilities  are  driven  by  two  mechanisms:  (i) 

resistive  phase  lags  due  to  magnetic  diffusion  and/or  (ii)  magnetic  repulsion 

between  the  beam  current  and  the  plasma  return  current.  The  first  mechanism  is 

in  some  ways  more  delicate  and  is  particularly  sensitive  to  the  details  of  a 

theoretical  model,  since  the  destabilizing  effect  of  resistive  phase  lags  between 

the  perturbed  beam  current  and  the  perturbed  currents  induced  in  the  plasma  can 

be  neutralized  by  the  stabilizing  effect  of  phase  mixing  among  beam  particles 

with  different  betatron  frequency. The  distribution  of  betatron 

frequencies  depends  on  the  degree  of  anharmonicity  of  the  self-pinch,  and 

therefore  on  the  radial  profile  of  net  current  density  .^(r);  all  particles  have 

the  same  betatron  frequency  if  Jno^r)  is  flat.  Thus  it  is  important  in  treating 

resistive  phase  lag  effects  to  model  the  beam  current  profile  J,  (r)  and  the 

bo 

return  current  profile  J  (r)  with  appropriate  rounded  shapes  and  to  include 

po 

phase  mixing  in  a  realistic  way.  When  this  is  done,  the  hose  mode  is 

found  to  be  unstable  even  if  there  is  no  equilibrium  return  current.  However  it 
is  widely  believed  that  for  typical  rounded  beam  profiles  such  as  the  Bennett 
profile  none  of  the  other  modes  are  destabilized  by  resistive  phase  lag  effects 
alone  when  phase  mixing  is  modeled  realistically;  these  modes  become  unstable 
only  when  the  return  current  exceeds  a  threshold  value. 

The  purpose  of  the  present  paper  is  to  investigate  the  destabilizing  effect 
of  plasma  return  current  on  a  variety  of  beam  modes  with  azimuthal  mode 
number  0  <  n  <  3  and  radial  mode  number  in  the  range  0  <  n  <  3.  We  exclude  the 

Manuscript  submitted  December  29, 1981. 
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effect  of  magnetic  phase  lags  (and  considerably  simplify  the  analysis)  by 

considering  only  frequencies  in  the  range  «  1,  where  is  the  dipole 

magnetic  decay  time.  Instabilities  in  this  frequency  range,  due  entirely  to 

plasma  return  current,  are  not  stabilized  by  phase  mixing,  and  are  only 

moderately  sensitive  to  the  profile  of  ( r ) ,  as  is  shown  in  Ref.  8  for  the 

sausage  mode  (m  =  0,  n  =  1).  We  consider  here  the  simplest  beam  equilibrium 

model,  in  which  J,  (r)  is  flat.  We  also  assume  here  that  the  conductivitv  of  the 
bo 

background  plasma  is  high  enough  to  assure  space  charge  neutrality.  Within  these 
constraints,  we  are  able  to  formulate  a  Vlasov-tlaxwell  treatment  of  the 
instabilities,  and  solve  it  in  closed  form. 

The  present  paper  is  a  direct  extension  of  our  previous  work^  to  a  wider 
selection  of  modes.  Analysis  of  the  higher  modes  considered  here  is  very 
complicated,  but  the  present  analysis  is  tremendously  simplified  by  the 
assumptions  of  charge  neutrality  and  «  1;  otherwise,  the  equilibrium 
considered  is  similar  to  that  of  Ref.  5.  With  slight  modifications,  the  results 
of  this  paper  can  also  be  applied  to  ion  beam  transport  in  an  inertial 
confinement  fusion  reactor  vessel. 

The  outline  of  the  paper  is  as  follows.  The  equilibrium  is  specified  in 
Sec.  II.  A  normal  mode  treatment  of  the  linearized  Vlasov-!laxwell  stability 
problem  is  formulated  in  Sec.  Ill  and  the  method  of  exact  solution  is  outlined. 
The  dispersion  relations  for  all  unstable  modes  with  0  <  u  <  3,  1<  m  +  n  <  3  are 
calculated  in  Sec.  IV  in  terms  of  the  degree  of  current  neutralization,  and  are 
discussed  with  emphasis  on  instability  thresholds  and  transitions  between 
oscillatory  and  purely-growing  instability. 
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II.  Equilibrium  Model 

We  consider  a  self-pinched  relativistic  electron  beam  propagating  through  a 
collision-dominated  background  plasma.  Cylindrical  polar  coordinates 
(r,0,z)  are  used,  with  the  z  axis  along  the  axis  of  symmetry.  Both  the  beam  and 
plasma,  in  equilibrium,  are  taken  to  be  azlmuthally  symmetric  (9/36  -  0), 
infinitely  long,  and  axially  symmetric  (3/3z  ■  0).  There  is  no  external  magnetic 
field,  and  we  assume  that  the  plasma  conductivity  o(r)  is  large  enough  to  assure 
space  charge  neutrality  out  to  a  radius  R^  much  larger  than  the  beam  radius  R^. 
This  requires  only  a  weak  conductivity, 


4iro(r)R.  >  c,  0  <  r  <  R.,  (la) 

D  D 

4*o( r)r  >  c,  R,  <  r  <  R  ,  (lb) 

b  c 

and  is  usually  satisfied  when  a  beam  propagates  into  pre-lonizcd  plasma,  or  even 
into  neutral  gas  of  density  >  1  torr  (except  at  the  beam  head).  We  also  assume 
that 


v/Yb  » 


V 


tic3(irb2-i)l/2 


17  (Yb-l)1/2kAnp 


«  1. 


(2) 


where  v  is  Budker's  parameter,  1^  the  bean  current,  -e  and  M  the  charge  and  mass 
of  the  electron,  c  Che  speed  of  light  in  vacuo,  and  Yb  the  bean  relativistic 
factor.  A  beam  satisfying  Eq.  (2)  is  paraxial,  i.e.  perpendicular 
velocities  v^  are  much  smaller  than  axial  velocities. 


We  consider  the  particular  equilibrium  distribution  function 


(3) 


£bo(H’V  -  2 7TS  *<"  -  ?  *2>  ‘<V  7  *V> 

D 

j. 

2  4  2  2  2 

where  H  *  (M  c  +  c  p  )  Is  the  particle  energy,  p  ■  p  -  eA  (r)/c  is  the  axial 

Z  Z  O 

canonical  momentum,  AQ(r)  is  the  axial  component  of  vector  potential  from  which 
the  equilibrium  self-pinch  field  B^r)  is  derived,  £  *  (Pr.Pg,Pz)  is  the 

q  222  1/2  ~ 

mechanical  momentum,  Y.  is  defined7  by  Y.  =  (1  +  P  /Me),  and  n.  and  y  are 

d  b  z  bo 

constants.  For  a  paraxial  beam,  Y  is  only  slightly  larger  than  Yfe*  As  a 
consequence  of  (2),  eAQ/c  «  pz  and  the  axial  velocity  vz  is  nearly  constant  and 
the  same  for  all  particles, 


v 

z 


v- 


(4) 


We  further  assume  that  the  equilibrium  plasma  return  current  J  (r)  has  the  same 

po 

profile  as  the  beam  current  J,  (r), 

DO 


J 

po 


(r) 


-  f  Jbo(r)’ 


(5) 


where  the  fractional  current  neutralization  f  is  a  constant  with  0  <  f  <  1.  From 
(2),  (3)  and  (5)  it  follows  that  the  equilibrium  beam  density  profile  is 


where 


nbo(r) 


nbo,  0  <  r  <  H 
0,  ^  <  r, 


(6) 


2  2  ~  —2 

V  “ 2c  -  V  Yb  “a 


(7) 


is  the  beam  radius  squared, 
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(8) 


“62  *KbV  °  '  f)  *  2'"b„'V  <■  -  f)/v 

is  the  betatron  frequency  squared,  and  ui  .  is  the  plasma  frequency. 

pb 

The  conductivity  profile  o(r)  docs  not  appear  explicitly  in  either  the 
equilibrium  or  perturbation  calculations  in  this  paper,  but  it  can  be  inferred 
from  conditions  (1)  and  (5)  that  o(r)  is  relatively  larpe  and  nearly  independent 
of  r  for  0  <  r  <  R^,  and  then  falls  to  a  much  smaller  value  consistent  with  (1) 
for  <  r  <  Rc« 

The  properties  of  this  equilibrium  and  related  equilibria  are  discussed  in 
more  detail  in  Refs.  4  and  5. 
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III.  Linearized  Vlasov  Analysis 

As  discussed  in  Sec.  I,  our  aim  in  this  paper  is  to  determine  the  stability 
properties  of  a  number  of  modes,  all  of  which  (except  for  the  hose  mode)  are 
destabilized  primarily  by  the  presence  of  an  equilibrium  return  current.  To 
eliminate  the  additional  destabilizing  (and  complicating)  effects  of  resistive 
phase  lags  due  to  magnetic  diffusion,  we  consider  only  nodes  that  satisfy 

I co|  T  =  xo(0)Rfe2  |m|/2c2  «  l,  (9) 

where  u  is  the  eigenfrequency  and  T.  is  "dipole"  magnetic  decay  time,  which 
actually  is  a  roughly  characteristic  decay  time  for  magnetic  field  patterns 
corresponding  to  any  of  the  first  few  perturbation  modes.  We  adopt  a  normal  mode 
approach  in  which  any  perturbed  quantity  $j(x  ,t)  with  azimuthal  mode  number  in  is 
represented  as 

♦,(x,t)  ■  #(r)exp{i(mO+kz-!ot)  J ,  (10) 

where  k  is  the  axial  wavenumber.  Introducing  the  variable'^  t  3  t  -  z/P.c,  this 

b 

Fourier  representation  is  rewritten  as 

♦^(x,  t)  -  *(r)exp(i(m6-uT-nz/$bc)] ,  (11) 

where  fl  =  u  -  kf^c  is  the  shifted  frequency  "seen"  by  a  beam  particle  (but  not  in 
the  sense  of  a  Loren tz  transformation). 

Conditions  (1)  and  (9)  together  Indicate  that  c/tv  is  long  compared  to  the 
beam  radius 

I «1  «  c/Rb,  (12a) 
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and  the  subsequent  analysis  clearly  indicates  that  when  (2)  is  satisfied 

|fl|  <  few  •  upbBb  «  c/Rfe  (12b) 

for  all  modes.  It  follows  from  (1)  and  (12a)  that  for  r  <  1^, 

4*o  »  |w|,  (13) 


which  guarantees  that  the  perturbed  beam  space  charge  field  is  also  completely 
neutralized  by  the  plasma.  In  the  context  of  (1),  (2),  (12),  and  (13),  only 
magnetic  fields  in  the  (r,0)  plane  play  any  significant  role  in  the 

A 

electrodynamics.  These  fields  can  be  derived  from  the  axial  component  A  of  the 
perturbed  vector  potential,  and  Maxwell's  equations  reduce  to  Ampere's  law 

A 

for  A.  Moreover,  (9)  indicates  that  the  perturbed  plasma  current  can  be 
neglected,  and  Ampere's  law  thus  reduces  to 


,  1  d_  d_ 
r  dr  r  dr 


2 

%  )  A(r) 


r 


4* 


Jb(r) 


(14) 


where  J. (r)  is  the  axial  component  of  the  perturbed  beam  current  density, 
b 

* 

Using  the  Fourier  representation  (11),  the  perturbed  beam  current  J  (r)  is 

D 

expressed  as  a  moment  of  the  perturbed  beam  distribution  function 
^bl^~l’®'*  T,z^» 


Jfe(r)  -  -  e  exp[i(-mfl  +  wr  +  «z/Bbc)J  /  d3pvzfbl(xit£,T,z),  (15) 


and  f  is  formally  evaluated,  in  standard  Vlasov  fashion,  as  an  Integral  over 
bl 


unperturbed  electron  trajectories. 


(16) 


fbl<&l’£  •T*z) 


_n  ,  ~  V  X  B_ 

rO  _  (  iilz  ~  ~1 

*  —  00 


r  lslz  i 
exp  l  -  J-  ) 

b 


In  Eq.  (16),  the  unperturbed  particle  trajectory  is  specified  as  a  function  of  z, 
with  =  x^(z-z),  v  =  v(z-z)  being  the  coordinates  in  the 

(r,9)  plane  and  velocity,  at  an  earlier  time,  of  a  particle  with  coordinates  and 

A  A 

velocity  x^,v  at  z.  It  is  assumed  that  Im£2  >  0.  We  note  that  B  =  V  x  A  ,  and 
use  (2)  and  (4),  from  which  it  follows  that  v,  «  v  ,  and  that  v  and  t  are 

i  2  Z 

essentially  constant  over  a  particle  trajectory.  We  can  then  reduce  Eqs.  (14)  - 
(16)  to  the  form 


“-f)  t737'3V-  V‘(t> 

r 


=  -  (2/Rb)  6(r  -  Rfa)  (  A  (Rb)  +  £1  I  (Rb,  Pq)  ] 


-  4y  2  M2  £1  u  2  0  (R  -  r)  (3I/3P2)  _ 

b  p  b  1  pi  Po’ 


(17a) 


where  I(r,p^)  is  the  orbit  integral 


I(r,p1)  =  i  /  qW( d <J>/ 2 tt)  J_^dz  A(r  )  exp[im(8  -  0)  -  i£lz/Pbc], 


(17b) 


and 


©  (x) 


1,  x  >  0 
0,  x  <  0 


(18) 


is  the  step  function.  We  have  introduced  polar  coordinates  (p^,$)  for  momentum. 
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and  have  defined 


po2(r>  5  v ,2"e2(Rb2  - r2)-  (19) 

When  supplemented  by  the  boundary  conditions  required  at  r  =  0  for  cylindrical 
symmetry. 


[dA/dr]r=0  =0,  m  =  0 


(20a) 


A(0)  =0,  m  >  1 


(20b) 


and  metallic  wall  boundary  conditions  at  radius  >  R^, 


A(Rc)  «  0, 


(20c) 


Eqs.  (17)  completely  specify  the  eigenvalue  problem  for  frequency  0. 

For  the  equilibrium  under  consideration,  with  J,  and  J  uniform  out 

bo  po 

to  r  =  R^,  the  unperturbed  electron  orbits  in  the  (r  ,  0)  plane  are  ellipses,  and 
can  be  specified  in  the  closed  form 


x  (  z  ) 


V'“b 


cos 


$  sin 


(2’a) 


'  ~  pi  “az  wsz 

y  (z)  -  v~77rr  sin  9  sin  (rt)  +  r  sin  9  cos  (« r)*  (21b) 

Yb  0  V  V 

Civen  the  simple  orbits  (21),  exact  solutions  to  Eqs.  (17)  and  (20)  are  found  in 
the  form  of  polynomials 
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Mr)  -  r"  1  ./W 23  . 
J-0  3  ' 


0  <  r  <  R^  (22a) 


A(r )  =  [tn(r/Rc)/in(Rb/Rc>]  f  &y 


Rb  <  r  <  Rc,  m  -  0, 

(22b) 


A(r)  -  Rb2n,(Rb2m  -  Rc2m)'1(rra  -  R^r"™)  l  Rfe  <  r  <  R£ ,  m  >  1. 

(22c) 

Hie  integer  n  is  the  radial  mode  number,  and  can  take  any  value.  The 
coefficients  a^  are  functions  of  ft  as  well  .as  of  the  mode  number  (n,n).  The 
determination  of  the  coefficients  aj  and  of  the  dispersion  relation  is  a 
straightforward  but  lengthy  algebraic  exercise,  for  which  the  degree  of 
complexity  increases  rapidly  with  increasing  mode  numbers. 

As  an  example,  we  shall  outline  the  solution  for  the  case  of  the 
axisymmetric  modes  with  m  »  0,  n  <  3.  The  eigenfunctions  are  of  the  form 


A(r) 


V 


2j 


0  <  r  <  Rb, 


3 

A(r)  -  [An(r/R  )/*n(R./R  )]  l  a 
C  *  C  >0  j 


Rb  <  r  <  Rc* 


(23a) 


(23b) 


We  first  substitute  (23a)  and  (21)  into  (17b)  and  evaluate  the  orbit 
integral  I  in  terms  of  ft  and  the  coefficients  a^.  This  part  of  the  calculation 
is  tedious,  but  can  be  carried  out  in  closed  form.  Next  we  substitute  (23)  into 
(17a)  and  solve  (17)  inside  the  beam.  Finally,  we  apply  an  appropriate  boundary 
condition  at  r  -  R^,  determined  by  multiplying  Eq.  (17a)  by  r  and  integrating 
over  R^  -  e  <  r  <  +  e,  with  e  ♦  0.  The  result  is  a  matrix  equation  of  the 

form 
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with  matrix 


where  the 


elements 


Xoo  =  Un(Rb/Rc)fl, 

*01  •  *00  +  2  ai  '  °» 

^2  =  %)  +  4C1  <6C2  +  °  '  4* 

Xq3  "  Xoo  +  651  <l20«253  +  12C2  +  l)  "  6* 

*11  "  *1  "  l’ 

X12  “  "  24^^, 

X22  B  ?!  <18^2  +  *>  ’  !» 

x23  “  *9/5)  (x22  ~  x33* 

x33  ”  V1200^  +  m2  +  °  "  1* 


(24) 

(25a) 

(25b) 

(25c) 

(25d) 

(25e) 

(25f ) 

(25g) 

(25h) 

(25i) 

(25J) 


's  are  susceptibility  functions,  given  by 
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{i  ■  -  vV(q2  -  * 


-  Ug2(fl2  -  16“g2 )_1 • 


c3  -  Wp^O2  - 


(26a) 


(26b) 


(26c) 


According  to  Eq.  (24),  the  dispersion  relations  for  the  modes  with  radial 
mode  number  j  are  specified  by  the  conditions 


0,  J  -  1,  2,  3. 


Since  Xjj  is  generally  a  polynomial  in  8  ,  modes  occur  in  pairs  +  8,  and  in 
general  there  are  several  distinct  pairs  of  modes  with  azimuthal  and  radial  mode 
numbers  (m, j).  Because  of  condition  (9),  we  have  effectively  specialized  to  the 
case  w  -  0;  thus  our  "dispersion  relation"  specifies  8  as  a  function  of  the  only 
other  free  parameter,  f.  To  find  the  eigenfunction  for  any  mode,  the  matrix 
equation  (24)  is  solved  for  the  a^,  using  the  previously  determined  value  of  8  to 
evaluate  the  matrix  elements 


IS 


IV.  Stability  Analysis 


A.  Results 


From  Eqs.  (24)-(27),  the  dispersion  relations  for  the  axisymnetric  (m 


0)  modes  are 


n  2  =  i  -  2f 


for  the  m  ■  0,  n  *  1  sausage  node. 


fl  4  -  (lOf  -  1)  fl  2  +  f  (16  f  +  1)  --  0 
e  e  e 


for  the  m  *  0,  n  »  2  hollowing  node,  and 


(fl  2  -  2f  )  (fl  2  -  8  f  )  (ft  2  -  18  f  )  +  (ft  2  -  8  f  )  (ft  2  -  18  f  ) 
e  e  e  e  e 


+  24  f  (fl  2  -  18  f  )  +  300  f  2  -  0 
e  e  g 


for  the  ra  *  0,  n  =  3  mode,  where  we  have  introduced  the  notation 


fe  i  1  -  £, 


■a  =  a/»  (j,  *  20fc/v 


In  similar  fashion,  we  have  calculated  the  dispersion  relations  for  the 
modes  (ra*l;  n*0,l,2),  (m«2;  n"0,l)  and  (m**3,n“0).  For  convenience,  we  list  all 


of  these  results  here: 


(1,  0):  fl  2  -  -f 


(hose  mode). 


IS 


(1,  1):  n  4  -  <5fe  -  1)  n  2  +  (9/4)fe2  -  0,  (34) 

(1,  2):  (2  n  2  -  9  fe)  (2  fl  2  -  25  f  )  +  4  «  2  -  0,  (35) 

(2,  0):  2  Q  2  =  4f  -  1,  (36) 

(2,  1):  2  »  8f  -  1,  (37) 

e 

(3,  0):  4  a  4  -  2( lOf  -  1)  n  2  +  3  f  (1  +  3f  )  =  0.  (38) 

e  e  e 


The  dispersion  relations,  i.e.  the  growth  rate  and  oscillation 
frequency  0^  as  functions  of  f,  are  plotted  in  Fig.  1  for  the  ra  =  0  nodes,  in 
Fig.  2  for  the  m  =  1  modes,  and  in  Fig.  3  for  the  r.i  =  2  and  n  ■  3  nodes,  and  uil 
be  discussed  in  the  next  section. 
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Fig.  1  —  Plots  of  normalized  (a)  growth  rate  £2j  and  (b)  shifted 
real  frequency  nr  versus  f  (Eqs.  (28)  —  (30)]  for  the  azimuthally 
symmetric  perturbations  (m  =  0). 


Fig.  2  -  Plots  of  normalized  (a)  growth  rate  S2j  and  (l>)  shifted  real 
frequency  nr  versus  f  (Eqs.  (33)  -  (35)1  for  the  m  =  1  modes  with 
radial  mode  number  n  =  0  to  2. 


Fig.  3  —  Plot*  of  normalised  (a)  growth  rate  ftj  and  (b)  shifted  real 
frequency  Slt  versus  f  [Eqs.  (36)  —  (38)]  for  the  (m,  n)  *  (2,  0), 
(2, 1)  and  (3,  0)  modes. 


B.  Discussion 


As  shown  in  Figs.  1-3,  for  any  given  mode  there  is  a  threshold 

value  f  such  that  instability  occurs  for  f  >  f^.  In  addition,  each  mode  has  a 

second  threshold  f^  >  f^  such  that  the  instability  is  purely  growing,  i.e. 

Reft  =  0,  if  f  >  f^.  If  f^  =  fj,  the  instability  is  always  purely  growing,  but 

if  >  fj  the  instability  is  oscillatory  in  the  range  f  <  f  <  A  purely 

growing  instability  mode,  in  the  limit  wt,  <<  1  which  we  have  considered,  may  be 

d 

expected  to  be  particularly  dangerous  to  beam  integrity.  The  values  of 
f^  and  f^  are  tabulated,  for  each  of  the  modes  in  Table  1. 

Among  the  ra  =  0  modes  the  n  =  2  hollowing  mode  would  appear  to  be 
particularly  dangerous  because  of  its  low  instability  threshold;  however  the 
instability  is  oscillatory  except  at  a  very  high  return  current  fraction, 
f^=  0.95.  The  n  =  1  sausage  node  has  a  somewhat  higher  instability  threshold, 
f^  =  0.5,  but  is  purely  growing  whenever  it  is  unstable,  and  has  the  largest 
growth  rate  of  any  mode  over  the  range  0.  6  <  f  <  1.  There  are  two  distinct  n  =  3 
branches,  one  of  which  is  purely  growing,  but  witli  the  large  instability 
threshold  f^  =  0.95,  while  the  other  is  always  oscillatory. 

Among  the  m  *  1  modes,  the  n  =  0  hose  mode  (which  has  been  studied  more 
extensively  in  Refs.  3-7)  would  appear  to  be  by  far  the  most  dangerous.  It  is 
unstable  for  any  value  of  f  >  0,  is  always  purely  growing  (when  «  1  and  thus 
resistive  phase  lag  effects  are  neglected,  as  is  done  throughout  this  paper),  and 
has  the  largest  growth  rate  of  any  m  *  1  mode  except  at  very  large  values  of  f. 
The  n  *  1  and  n  ■  2  modes  are  oscillatory  up  to  large  values  of  f.  Both  of  these 
modes  turn  on  at  fj  ■  0.5,  and  have  the  same  growth  rate  in  the  regime  where  both 
are  oscillatory. 

The  m  ■  2  modes  are  elliptical  distortions  which  at  large  amplitude  tend 
to  split  the  bean  in  two.  Both  the  (n,n)  ■  (2,0)  and  (2,1)  modes  are  purely 


growing  whenever  they  are  unstable,  but  the  threshold  for  both  modes  is  rather 
high,  fj  *  0.75  and  0.87  respectively.  The  (m,n)  =  (3,0)  f ilaraentation  mode,  on 
the  other  hand,  turns  on  at  the  surprisingly  low  threshold  f}  =  0.53,  but  is 
purely  growing  only  at  very  large  f. 

Of  all  the  modes  studied,  the  dispersion  properties  would  appear  to 
single  out  the  hose  as  most  dangerous,  with  two  of  the  axisymnetric  nodes,  m  **  0, 
n  =  1,  2,  next.  This  is  in  accord  with  the  conventional  wisdom  that  hose  is  the 
most  dangerous  mode.  On  the  other  hand,  the  hose  mode  in  a  well-prepared  beam 
must  grow  out  of  small-amplitude  noise,  whereas  the  axisymnetric  nodes  may  be 
initially  excited  at  much  larger  amplitude,  due  to  the  violent  pinch— down  of  the 
beam  head^  * or  to  the  axisymnetric  oscillations  of  a  beam  which  is  mismatched 
at  injection.  Thus  experimental  study  of  all  three  of  these  modes,  as  well  as 
further  theoretical  investigation  using  simulation  techniques  to  model  the  beam 
realistically,  would  appear  to  be  warranted. 
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Table  I 


Mode  Numbers 

m 

n 

*1 

{2 

0 

1 

0.5 

0.5 

0 

2 

0.38 

0.95 

0 

3 

0.45 

1 

0 

3 

0.95 

0.95 

1 

0 

0 

0 

1 

1 

0.5 

0.89 

1 

2 

0.5 

0.96 

2 

0 

0.75 

0.75 

2 

1 

0.87 

0.87 

3 

0 

0.53 

0. 96 

Threshold  values  of  fractional  current  neutralization  f  for  the  various 
modes:  f  >  f ^  Is  the  condition  for  instability,  and  f  >  f^  is  the  condition  for 
pure  growth. 


I 


C.  Concluding  Remarks 

In  this  paper  we  have  presented  a  stability  analysis  of  all  of  the  first 
nine  modes  (1  <  a  +  n  <  3)  of  a  beam  propagating  in  a  resistive  plasma.  The 
results  are  based  on  an  exact  Vlasov  analysis,  which  is  carried  out  in  closed 
form.  This  is  possible  because  of  two  principal  restrictions:  we  considered 
only  beam  equilibria  in  which  JbQ(r)  and  JpQ(r)  are  flat»  and  we  considered  only 
the  case  u»xd  «  1.  These  two  restrictions  are  compatible  with  each  other:  in 
the  regime  «  1  the  instability  is  driven  only  by  the  presence  of  an 
equilibrium  return  current  1^  *  0,  and  this  mechanism  is  only  moderately 
sensitive  to  the  profile  JbQ(r).  The  analysis  is  obviously  incomplete.  Modes 
with  ux^  ~  1  are  also  important,  and  because  these  modes  are  driven  additionally 
by  resistive  phase  lags  between  perturbed  beam  current  and  perturbed  plasma  eddy 


larger  growth  rates.  However  these  modes  are  also  sensitive  to  the  stabilizing 

effect  of  phase  mixing  among  particles  of  different  betatron  frequency,  and 

therefore  to  the  profiles  of  JbQ(r)  and  JpQ(r).  Thus,  if  one  is  to  treat  these 

instabilities  over  the  full  range  of  wt, ,  it  is  important  to  use  realistic  beam 

d 

profiles  and  not  to  use  models  that  eliminate  or  misrepresent  phase  mixing.  A 
complete  Vlasov  analysis  would  of  course  be  ideal,  but  no  such  analysis  has  been 
presented  to  date  for  beams  with  rounded  profiles.  Instead,  several 
phenomenological  techniques  have  been  developed  to  model  phase  mixing  and  beam 
profile  dependence  of  the  hose®’^*®^  and  sausage®*®  nodes,  but  these  techniques 
have  not  been  adapted  to  the  higher  modes  discussed  in  this  paper,  as  yet. 

One  other  restriction  of  the  present  model  should  be  mentioned.  We  have 
implicitly  assumed  that  the  conductivity  channel  o(r)  is  pre-formed  before  the 
arrival  of  the  beam  and  is  uninfluenced  by  the  passage  of  the  beam.  In  fact,  for 
a  high-current  beam  injected  into  a  neutral  or  weakly  ionized  gas,  augmentation 
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of  the  conductivity  as  a  result  of  ionization  and  heating  by  the  beam  can  have  a 
strong  stabilizing  effect, ^  and  should  be  treated  self -‘consistently  with  the  beam 
dynamics  in  these  situations. 
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